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Abstract 

The Poincare Equivalence Theorem states that any optical element which contains no absorbing 
components can be replaced by an equivalent optical model which consists of one linear retarder and 
one rotator only, both of which are uniquely determined. This has many useful applications in the 
field of Optics of Polarized Light. In particular, it arises naturally in attempts to reconstruct spatially 
varying refractive tensors or dielectric tensors from measurements of the change of state of polarization 
of light beams passing through the medium, a field which is known as Tensor Tomography. A special 
case is Photoelasticity, where the internal stress of a transparent material may be reconstructed 
from knowledge of the local optical tensors by using the stress-optical laws. - We present a rigorous 
approach to the Poincare Equivalence Theorem by explicitly proving a matrix decomposition theorem, 
from which the Poincare Equivalence Theorem follows as a corollary. To make the paper self-contained 
we supplement a brief account of the Jones matrix formalism, at least as far as linear retarders 
and rotators are concerned. We point out the connection between the parameters of the Poincare- 
equivalent model to previously introduced notions of the Characteristic Parameters of an optical 
model in the engineering literature. Finally, we briefly illustrate how characteristic parameters and 
Poincare-equivalent models naturally arise in Photoelasticity. 

Keywords: Poincare Equivalence Theorem, Matrix Decomposition, Characteristic Parameters, 
Equivalent Optical Models, Jones matrix formalism, Photoelasticity 



1 Introduction 

The propagation of an electromagnetic wave through a material medium is the result of the interaction 
between the fundamental fields E and B with a macroscopic number of microscopic sources constituting 
the bulk matter. In principle, these sources must be incorporated into the dynamics of the total system 
by appropriate interaction terms. However, if the frequency bandwidth of the light under consideration 
is sufficiently far away from the resonance frequencies of the macroscopic medium, any photon-atom 
encounter is only transient; in this case we may refer to the medium as a passive optical element. 
Its macroscopic effect on impinging radiation may be summarized by introducing a (space- and time- 
dependent) refractive index. 

There are essentially two classes of such optical elements: Firstly, polarizers, which absorb, or at least 
attenuate, one of two given orthogonal polarization forms of the light passing through the medium; and 
secondly, retarders, which introduce a phase retardation between the components of these polarization 
forms, but otherwise preserve the total intensity of the light. In principle, a polarized light beam can 
be analyzed with respect to any two orthogonal polarization forms. However, in practice, two cases are 
of particular importance: The linear retarder, which introduces a phase lag between the components of 
linearly analyzed light, and the rotator, which does the same with respect to the circularly polarized 
components of a polarized light beam. 
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The description of the passage of light through optical elements depends on whether the light is polar- 
ized or not. For unpolarized or partially polarized light, a description in terms of the Stokes parameters 
[1] (for a definition, see e.g. [2]) is appropriate; if these are arranged into four-dimensional real vectors 
such that optical elements act on it by means of real (4, 4) matrices we speak of the Mueller formalism 
[3, 4]. If, on the other hand, the light beam is completely polarized, the polarization state of the trans- 
verse degrees of freedom of the electric field can be conveniently arranged into a two-dimensional complex 
vector, while optical elements now act on it by means of complex matrices. This formalism is called the 
Jones calculus [5, 6, 7, 8, 9, 10, 11, 12, 13]. The above-mentioned concepts prove to be useful not only 
within a classical context, but appropriate generalizations are used in Quantum Optics as well: The 
unitary representation of beam splitters in the Jones matrix formalism acting on quantum-mechanical 
mode operators was discussed in [14]. A broad introduction to optical elements in linear optical networks 
within a quantum-optical context was recently given in [16]. 

Within the framework of the Jones calculus, polarizers are described by Hermitean matrices, while 
retarders are represented by unitary ones. The latter reflects the fact that retarders preserve the total 
intensity of the light. The fact that unitary matrices have the group property then implies that an 
arbitrary sequence of retarders is again described by a unitary matrix. Such a sequence could be a pile of 
bircfringent plates, each with constant rotation of principal axes and constant phase retardation. But we 
can also consider the limit of infinitely many, infinitely thin, bircfringent plates: This situation describes 
a medium which exhibits a spatially varying dielectric tensor , or more generally, refractive tensor n,j . 
In this case, the phase velocity depends not only on the location within the medium, but also on the 
polarization state of the (locally-plane) wave at this point. This corresponds to the most general case of 
an inhomogeneous anisotropic medium. 

Such a scenario arises naturally in the field of Photoelasticity [17, 18, 19, 20, 21, 22, 23, 24]: Experience 
shows that certain materials, such as glasses and polymers, are optically isotropic and homogeneous when 
unloaded, but exhibit local anisotropy when strained by an external load. The relation between the 
resulting dielectric tensor and the stress tensor cry in the interior of the medium is called the stress- 
optical law; its basic form has been discovered long ago by Maxwell [25]. Integrated Photoelasticity is 
concerned with the reconstruction of optical tensors, and via the stress-optical law, of stress tensors, from 
data sets which are acquired by sending polarized light through the loaded specimen at many different 
angles, thereby measuring the change in the state of polarization. For this method to work, the material 
must be sufficiently transparent to ignore absorption within the medium. Since tensorial quantities are 
reconstructed, these methods belong to the field of Tensor Tomography. 

In the examples above, the medium is given by a pile of retarders and/or rotators, possibly in the 
infinitesimal limit. The previous discussion implies that, within the Jones formalism, such a medium may 
be described by a unitary (unimodular) matrix U. It turns out that such a matrix always determines 
- and, in turn, is fully determined by - a set of characteristic directions [26] and a characteristic phase 
retardation between them. These characteristic parameters [26] have been introduced within the context 
of Photoelasticity from an engineering point of view, and are not necessarily standard in the mathematical 
literature [a rigorous mathematical definition will be given in section 3]. If the matrix U is interpreted 
as belonging to an optical element of the kind as discussed above, the characteristic parameters can be 
given an operational meaning: The primary characteristic directions determine those planes of linear 
polarization at the entry into the medium for which the state of polarization at the emergence from the 
medium is again linear. The secondary characteristic directions determine the planes of linear polarization 
of the emerging light, if the incident light was linearly polarized in the primary characteristic directions; in 
general, they differ from the primary ones. It turns out [sections 3, 7] that there are always two orthogonal 
primary and two orthogonal secondary characteristic directions. Light which is linearly polarized along 
the two primary directions travels with different phase velocities, though, so that both waves emerge with 
a phase difference - the characteristic phase retardation. 

It is an important task to reconstruct the optical properties of the medium, i.e. the matrix U , from 
measurements of the characteristic parameters. To this end, relations between the parameters of a normal 
form of unitary matrices and the characteristic parameters may be derived [22, 23, 26]. However, these 
relations only determine the squares of sines or cosines of angles and therefore the actual computation of 
these quantities is rather involved. It is at this point where the Poincare Equivalence Theorem [27] comes 
into play: This theorem states that any optical medium which is described by a unitary unimodular 
matrix can be replaced by an optically equivalent model which consists of one retarder followed by 
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one rotator, or the other way round; in each case, both elements are uniquely determined. Then it 
can be shown easily that the characteristic parameters of the optical medium coincide with the optical 
parameters of the equivalent model, i.e. with the principal directions and phase retardation for the 
equivalent retarder, and the rotation angle for the equivalent rotator; but from the latter, the matrix U 
may readily be determined completely. This reconstruction becomes particularly important in the case 
of tensor tomography of photoelastic media, where the overall effect of the specimen on the passage of 
light represents the accumulated effect of many (infinitesimal) retarders and rotators. 

These examples highlight the importance of the Poincare Equivalence Theorem in the field of Optics 
of Polarized Light. Although references to the equivalence theorem, or applications thereof, are made 
frequently in the literature [20, 22, 23, 24, 28, 29], its precise content is apparently not widely published; 
a rigorous proof of the theorem, or some equivalent statement thereof, seems even less readily available. 
This observation constitutes the starting point of our paper: We rigorously prove a matrix decomposition 
theorem for unitary unimodular (2, 2) matrices [section 4], from which the Poincare Equivalence Theorem 
follows as a corollary [section 5]. To put this result into a context, and to make the paper reasonably 
self-contained, we provide a brief overview of the description of polarized light, and optical elements 
acting on it, in terms of the Jones calculus; this is done in section 2. The second objective of this paper 
is to clarify the relation between the (operationally defined) characteristic parameters of a transparent 
model according to Aben [22, 23, 26], and the optical parameters of the Poincare-equivalent model. 
To this end we recapitulate in section 3 the notions of characteristic directions and characteristic phase 
retardation for non-absorbing optical elements. We also prove the existence of characteristic directions for 
any unitary unimodular (2, 2) matrix in a standard parametrisation. The relation between characteristic 
parameters and the parameters of the Poincare-equivalent retarder and rotator then are analyzed in 
section 7. Finally, in section 8 we outline the basic ideas of three-dimensional photoelasticity in order to 
illustrate the context in which a specimen with continuous variation of optical tensors presents itself, and 
how the Poincare-equivalent model may be important to the objective of tensor tomography. In section 
9 we summarize our results. 



2 Polarized Light in the Jones Matrix Formalism 

In this paper we only consider completely polarized classical light. Strictly speaking, this is an idealization, 
which can be approximately realized only if 1) the light is strictly monochromatic, and 2) the complex 
degree of coherence [2] of two orthogonal components of the electric field always has unit modulus. 
Physically, this means that the two orthogonal components have a well-defined and constant phase relation 
for all times. 

The transverse electric field of a general eliptically polarized light beam propagating in the 3-direction 
takes the form 

E(x, t) = a\ cos {kz — ut + 5\) ei + a 2 cos (kz — ut + 5 2 ) e 2 , (1) 

where a\,a 2 are real amplitudes, u = uj/k is the phase velocity in the (isotropic and homogeneous) 
medium, 5\ and 8 2 are constant phases, and ei and e 2 are real unit vectors in the direction of the x- and 
y-axes. Since the action of optical elements on light beams is accomplished by linear transformations, it 
is admissible to assume that the beam has unit relative intensity [21] a\ + a\ = 1. We now choose a fixed 
location within the beam, say, z = 0, and express (1) as the real part of the complex vector 

E = e~ iujt { oi e iSl ei +a 2 e 2 } . (2) 

The quantity in curly brackets defines the polarization form of the light; the associated complex 2-vector 

JSi" 



= rf«j . (3) 



v 2 j \a 2 e 

is called the Jones vector of the light beam with respect to linear polarization along the directions ei 
and e 2 . We use a subscript I to denote the fact that the Jones vector (3) denotes linearly polarized 
components. 

As can be seen, the polarization form and its associated Jones vector are defined only up to a global 
phase, i.e. Jones vectors v and e l ^v are equivalent, since the global phase $ can always be reabsorbed 
into the total phase e~ lu)t governing the time dynamics. 
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Eq. (2) represents the decomposition of the (complex form of the) electric field vector into two 
linearly polarized components whose amplitudes are ai,a%, and whose relative phase is 5 = 5% — 5\- 
This decomposition is convenient to derive the Jones matrix of a linear retarder. To study the rotator 
it is necessary to analyze the light with respect to complex basis vectors e + and e which represent 
right-handed and left-handed circular polarization, respectively, which are defined as 

e± = -j= (ei ± ie 2 ) , (4a) 
(e+ , e_) = (ei , e 2 ) M , M = _fj . (4b) 

[We recall that right-handed = left-circular, and vice versa]. The electric field vector (2) can be expanded 
in terms of the basis (e + ,e_), 

/y i \ - („ „ ^ ( v + 



E~( ei ,e 2 ) (^j=(e + ,e_) , (5) 

such that the components v + and w_ make up the Jones vector v c of the light in the basis (4), 

Dfi) • (6) 

We use a subscript c to point out that the Jones vector (6) denotes circularly polarized components. 

Evidently, the matrix M mapping the basis vectors (ei,e2) for linear polarization to those (e+,e_) 
for circular polarization is unitary. Indeed, this is true for any transformation between basis vectors of 
different polarization forms. 

We are now in a position to derive the Jones matrices for linear retarders and rotators. 



2.1 The linear retarder in the Jones calculus 

A linear retarder is an optical clement made of a material which exhibits homogeneous optical anisotropy 
along a given direction of light transmission [2, 20, 21]. Perpendicular to the direction of passage of 
light, the retarder has two distinct orthogonal directions, called the fast and the slow axis, respectively. 
The phase velocity of a plane wave which is linearly polarized along the fast axis is greater than for 
waves polarized along the slow axis. This means that, at the point of emergence from the retarder, 
the component along the slow axis ("slow component") has acquired a phase lag 5 with respect to the 
component along the fast axis ("fast component"). Let us assume that the fast /slow axis is oriented 
along the x/y-axis. If the Jones vector of light prior to entry into the retarder is 



V C2 , 

then after emergence from the retarder it has changed to 

The exponential e lS indeed represents a phase lag of the local harmonic oscillator at z = vibrating in the 
y-direction, because of our convention to count the phase of the time dynamics as e~ luJt . Furthermore, 
<fr g represents a global phase which is picked up by any wave passing through the retarder, irrespective 
of its polarization form. Since global phases arc immaterial in our present discussion, it is common to 
adjust them for maximum convenience, or otherwise discard them altogether, which is what we shall do 
throughout. Physically, a global change of phase could be accomplished by passage of the light through 
an isotropic medium whose optical path length corresponds to the desired phase difference. 

Since the global phase in (8) can be chosen freely, we may extract and discard a global phase <& s = |, 
in which case the transformed Jones vector reads 



e - 5 / 2 ci 



(9) 
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It is now clear that the matrix which accomplishes the transformation v — > v' must take the form 

s e -i6/2 o 



Jlin(0,J) 







JS/2 



(10) 



Our notation Ju n (0, 5) reflects the fact that the phase difference between the components is 5, while the 
angle between the fast axis and the x-direction is 6 = 0. 

It is easy to derive the Jones matrix Jii n (0, 5) of a linear retarder whose fast axis makes a nonvanishing 
angle with the £-axis. To this end we introduce the rotation matrix 



R(9) 



cos 9 sin 9 
— sin 9 cos 9 



(11) 



This is an SO(2)-matrix representing a passive rotation by 9 about the z-axis, that is to say, the basis 
vectors and vector components transform according to 



E ~ (ei , e 2 ) 



= K , e 2 ) 
= R(6) 



(ei, e 2 ) = (e 1; e 2 )R(9) T 
R(9) T = R(9)- 1 = R(-9) 



(12a) 
(12b) 



The group property implies that a succession of two rotations is again a rotation; furthermore, two 
rotations around the same axis commute. Thus we have 



R(0i + 2 ) = R(0i)R( 



R(0 2 )R(0i) 



(13) 



Using the rotation matrices it is easy to derive the Jones matrix Jij n (#,<5) for a linear retarder with 
nonvanishing angle 9: We only need to keep in mind that in the reference frame of the fast (ei) and slow 
(e 2 ) axes, the retarder must act with a Jones matrix (10) on the components v[ and v' 2 . Thus, 



E ~ (ei , e> 2 ) 



(ei , e 2 )J lin (0,<5) 



= (ei , e 2 )R(9) T 3 Via (0,S)R(9) 



= (ei , e 2 ) Jun(M) 



where we have used eqs. (12). It follows that 

Ji in (M) = R(-0)Jim(o,<5)R(0) 

or explicitly, 



/ cos i sin | cos 29 

JlinOM = 2 . . ] . „ fl 

V — % sin 2 sin2C 



COS ; 



-i sin | sin 29 N 
+ sin| cos 29 



(14) 



(15) 



(16) 



Equation (15) is indeed valid for Jones matrices J of arbitrary optical devices whose preferred axes 
are related by an orthogonal transformation, since all we have used in the derivation of (15) was the law 
of passive coordinate transformations onto a coordinate system whose axes coincided with the preferred 
axes of the optical device. 

We note that the Jones matrices of the linear retarder, (10) and (15), are unitary unimodular matrices, 
since we have chosen the global phase $ s in (10) in such a way that det Jn n (0, S) = 1. 



2.2 The rotator in the Jones calculus 

The rotator is an optical element which introduces a phase lag S of the left-handed relative to the right- 
handed component of the complex electric field vector (2). If the electric field is expanded in terms of 
the basis (4), the left-handed / right-handed components are just the elements i>_ and v+ in the Jones 
vector v c given in eq. (6). Thus, after passage through the rotator we have 

■ (17) 
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Evidently, the linear map which accomplishes the transformation (17) is precisely the matrix Jn n (0, 6) of 
the linear retarder given in (10); thus, in the basis (4) we simply have 



v' c = Jl in (0,(5) V c 



(18) 



Using eq. (6) which effects the transformation between the components for linear and circular polarization, 
we can easily transform eq. (18) back into the Cartesian frame (ei, e 2 ), in which the Jones vector is given 

by (3), 

w,' = MJi in (0,<y)Mtt; I =J Iot (-f)t; I . (19) 
This determines the Jones matrix J rot (— §) of the rotator in the Cartesian basis, 



*J rot ( 2 ' 



COS ; 



. sin ■ 




R( 



2> 



(20) 



where the rotation matrix R was defined in (11). 

The reason for our seemingly strange notation convention becomes clear if we recall that a positive 
phase lag of the left-handed component with respect to the right-handed component should rotate the 
plane of polarization of a linearly polarized light beam in a counter-clockwise sense; this is a positive 
angle +s^L in the xy-plane, if we look into the oncoming wave, i.e. from positive towards negative 
z-values. It is easy to check this: 



Jrot(-f) 



cos 9 
sin 9 



cos (9 + |) 
sin (9+1) 



(21) 



Here, [cos6*,sin#] T is the Jones vector of a normalized, linearly polarized beam, and 9 is the initial angle 
of inclination of the plane of polarization. After passage through the rotator, the plane of polarization 
has rotated by +6/2. 



This finishes our discussion of linear rctarders and rotators in the Jones calculus. 



3 Characteristic parameters of 577(2) matrices 

As mentioned in the Introduction, the terms 'Characteristic Directions', 'Characteristic Phase Difference' 
and 'Characteristic Parameters' were first introduced by Aben within the framework of Photoclasticity 
[26]. We now wish to explain their meaning. 

Consider a non-absorbing optical medium whose Jones matrix U is unitary. If U happens to have 
a determinant different from 1, we can always extract and discard a global phase, so that U may be 
assumed to be unimodular, det U = 1, and hence an element of SU(2). A parametrisation of SU(2) 
matrices is given by [30] 

fcos9e lcl> —s'm9e tx \ „ r _ _ n . r _ _ . ,_„. 

U =(^4 e -ix cosfle-*' ' 0G [°'i] ' &Xe[0,27r) • (22) 



sin 



This parametrisation is unique in the interior of the interval [0, -|], i.e. for < 9 < \. If 9 takes one 
of the boundary values, parts of the parameter space must undergo identification in order to make the 
map between parameters and matrices unique: For 9 = we must identify points (9, <j>, %) = (0,0, x), 
X € [0, 27r). For 6 = § , the points (9,cj),x) = (f ,<f>, 0), <f> € [0, 2n), must be identified. This agreement 
makes the mapping between parameter space and matrices bijective. 

Using this parametrisation we prove the 

Theorem 3.1 (Characteristic directions). For every unimodular unitary matrix U of the form (22) 
there exist two real vectors 

^=(1;) ' - =1 ' 2 ' ^ 
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such that the vectors U w m are real vectors times an overall phase, 

_ ( cos l'm\ _ J<S>„ 



Uw m = e 1 ^ \~£J£ j = e**~ w' m , m = 1, 2 . (24) 

77ie angles -f m are solutions of the equation 

. o sin(0 + x)sin26> 

tan2 7m = . . . , (25) 

cos ft sin 20 — sin ft sin 2% 

7ie[-f,f] a ^ 72=7i + f • (26) 
TTie vectors w m , m = 1, 2 as jtuen in fi?3) span £/ie primary characteristic directions o/ ifte matrix U . 
The angles j' m are solutions of the equation 

tan 2y = sin(0- X )sin2ft ^ 

cos 2 sin 20 + sin sin 2% 

7i€[-f,f] and 7 2 = 7i + f ■ (28) 
The vectors w' m ,m = 1,2 as gwen in span i/ie secondary characteristic directions of the matrix U. 
The angles 7 m and -f' m are well-defined if and only if 

(0,0) ^{ (0,0), (0,f), (0,tt), (0,f ) } , (29a) 
(6, X )?{ (f,0), (f,f), (|,7t), (f, f)} . (29b) 

J/ one of the parameter values (29) is taken, the right-hand side of eqs. (25, 27) becomes the ill-defined 
expression jj, and no angles j m , j' m can be determined. 

Proof: 

If the statement (24) holds, then the quotient w' ml /w' m2 , or equivalcntly, the product w' ml {w' m2 )* , 
must be real, for m = 1,2. This requires that 

3mw' ml ( W ' m2 )* =0 . (30) 

We now express («4i)i,2 as (Uw m )i^ 7 respectively, and impose the condition (30). Under the conditions 
(29) this leads to eq. (25), which determines an angle 271 lying between — ^ and Then, the angle 
271 + 7r also satisfies (25). It follows that 71 lies between — j and +j, and the second solution is 
72 = 7i + § • Thus, there are two primary characteristic directions, and they are orthogonal. 

To determine the angles Y m , we invert (24), 

U*v/ m = e- i *'»w m . (31) 

This has the form of eq. (24), and we can therefore use the same line of arguments as before, if we put 
the matrix [/t into the same form as (22), 

U^Ui => 0-»-0, 0^-0, X^X • (32) 

If the replacements (32) are used in (25) we obtain (27), again under the conditions (29). — It is easy 
to check that, whenever one of the values in (29) is taken, the right-hand sides in both (25) and (27) are 
ill-defined. — This finishes the proof. ■ 

The physical meaning of this is as follows: The existence of two solutions for (24) implies that there 
exist two perpendicular directions for the plane of linearly polarized light, given by the solutions (23), 
such that if light polarized in this direction enters the device it will emerge linearly polarized again; this 
is a consequence of the fact that the components of the transformed Jones vector have the same overall 
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phase. This statement is true only for the light as it emerges behind the apparatus; in the interior of 
the device, the light is elliptically polarized in general. Furthermore, the phase velocity of a linearly 
polarized wave entering the device at the angle 71 differs from the one entering at 72, so that both waves 
will emerge with a phase difference A, see eq. (36) below. 

We now turn to the determination of the phases <& m in (24) . We find 

Proposition 3.2 (Determination of phases). The phases $ m; m = 1,2 satisfy the equations 



and 



cos 2 $ m 


= cos 2 6 


cos 2 4> + sin 2 9 cos 2 x 7 


(33a) 


sin 2 $ m 


= cos 2 9 


sin 2 (j) + sin 2 9 sin 2 \ , 


(33b) 


tan 2 $ m 


sin 2 (f. 


1 + sin 2 x t an2 


(33c) 


cos 2 4 


< + cos 2 x tan 2 9 


cos2$ m 


= cos 2 9 


cos 2(f) + sin 2 cos 2x ■ 


(34) 



Proof: 

We perform the matrix multiplication Uw m with the form (22) for U, but the angles 7™ need not 
be known explicitly. On the right-hand side of this equation we have e 4 *" 1 w' m . We now expand this 
equation into real and imaginary parts, obtaining two equations which on their right-hand sides contain 
cos $ m w' and sin $ m w', respectively. If the component equations are squared and added, we obtain the 
relations (33a, 33b). To arrive at (33c), we divide (33a) by (33b) and then divide the whole equation by 
a common factor cos 2 9. Subtracting (33b) from (33a) then yields (34). ■ 

It is important to point out that the phase changes $ m in general do not coincide with the real phases 
of the light as it passes through the real device, for in the latter case there is also an additional global 
phase which is accumulated by any light beam that travels through the apparatus. This global phase 
has been discarded in the construction of the Jones matrices of the retarder and the rotator, in order 
to make these matrices unimodular. However, the global phase is the same for light of any polarization 
form, and therefore cancels out if we construct the phase difference between light beams which enter the 
device along the primary characteristic directions. This phase difference docs not depend on unobservable 
parameters and can therefore be measured. In principle, it can be evaluated from (33c): This equation 
has four solutions, given by 

$G [-§,§] , <S> + TT , -$ , -<P + TT . (35) 

One of the solutions $1 and $2 must be chosen from the set {$, $ + ir}, while the second one must take 
values in {—$,—$ + n}. In both cases we arrive at the conclusion that 

A = 2$ 2 mod 7T . (36) 

At this point a remark concerning the results (28) and (36) are in order: 

Eq. (27) determines the directions of polarization of the emerging light. However, it docs not deter- 
mine whether Uwi is equal to or w' 2 , and similarly for Uw2- This is reflected in the indeterminacy of 
the phase difference in (36): The modulus of A is determined only up to 7r, while the sign is undetermined. 
It is at this point where the decomposition theorem 4.1 and the ensuing Poincare Equivalence Theorem 
prove very useful, since they allow to remove this indeterminacy completely, provided the measurements 
on the sample are capable of distinguishing between the fast and the slow axis of the equivalent retarder 
[see sections 5, 7]. 

For every optical device described by a unitary unimodular matrix U, the set of quantities (71,7^ A), 
are called the characteristic parameters of the model. 

4 Matrix Equivalence and the Poincare Equivalence Theorem 

In theorem 4.1 we show that every two-dimensional unitary unimodular matrix U has a unique decom- 
position in terms of a minimal number of 'basic building blocks'. There are two of these building blocks: 
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The first is given by the rotation matrices R(£) as defined in (11); the second is given by diagonal unitary 
unimodular matrices 

V 5 



DW^O e - iS ) , ^e[0,27r) . (37) 
After that, in theorem 5.1 we deduce the Poincare Equivalence Theorem from theorem 4.1. 

Theorem 4.1 (Decomposition). Let U be a unitary unimodular matrix in the parametrisation (22). 
Then the following statements are true: 

I.) If the angle 9 in (22) lies in the interior of the interval [0, §], < 9 < f , then there exist angles 

a,/3e[0,7r) , (5e[0,27r) , (38) 

such that U has a decomposition 

I/ = R(-2±£)D(J)R(2=2) (39a) 



2 /"V"/"V 2 

sin^M 




'cos 6 cos /3 + i sin 5 cos a — cos 8 sin /3 + i sin <5 sin a 
cos (5 sin /? + i sin 5 sin a cos 5 cos /3 — i sin 5 cos a 



(39b) 
(39c) 



77ie angles (38) are uniquely determined, and there is no way to represent U by less than three 
matrices of the kind (11) and (37). 

II.) For 9 = we have 

U = D(0) . (40a) 
This is a unique minimal decomposition in terms of one diagonal matrix (37). 
A representation in terms of three matrices, form (39a), may be given for parameter values 

a =0 = , 5 = <j) , (40b) 

which is unique provided that 

<A£{0,f,7T,f } . (40c) 

If (p G {OjfjTT,^-}; o, decomposition in terms of three matrices is not unique, since one of the 
angles a or is left unspecified. However, the choice of parameters (40b) is unique in that it always 
corresponds to the minimal decomposition (40a). 

III.) For 9 = j we have 

f/ = R(-f)D(- X ) . (41a) 

This is a unique minimal decomposition in terms of one diagonal matrix ( 37) and one rotation 
matrix (11). 

A representation in terms of three matrices, form (39a), may be given for parameter values 

a = 0=l , S = -x , (41b) 

which is unique provided that 

X^{0,f,7T,f} . (41c) 

If X € {0,f,7T, a decomposition in terms of three matrices is not unique, since one of the 

angles a or is left unspecified. However, the choice of parameters (41b) is unique in that it always 
corresponds to the minimal decomposition (41a). 

Proof: 

We first must confirm that the right-hand sides of (39a, 40a, 41a) indeed lie in SU(2). But this is 
clear, since the adjoints of these matrices are equal to their inverses, and the determinant is equal to 
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one in each case. So, we know that these matrices comprise a subset of SU(2); but we do not yet know 
whether this subset comprises the whole of SU(2). This is what we are going to show next. 

To this end we first focus on decompositions in terms of three matrices, (39a). Let {6,4>,x) be a 
parameter triple pertaining to a given form (22). If a decomposition (39a) of U exists, there must exist a 
triple (a, (3, S) such that the real and imaginary parts of each component of the matrices (22) and (39c) 
coincide. This yields four equations 

cos S cos (3 — cos 9 cos <j> , (42a) 
cos 8 sin (3 — sin 9 cos x > (42b) 
sin 5 cos a = cos 9 sin <\> , (42c) 
sin S sin a = — sin 9 sin \ . (42d) 

We now show that the system (42) is indeed solvable for the triple (a, (3,6), and that the solution is 
unique under the circumstances described in case I.) - III.) above. 

We first examine case I.): 

On dividing (42c) by (42a), and (42d) by (42b), we find 

cos a 



tan </> = tan S 
tan x = — tan 5 



cos (3 
sin a 
sin/3 



(43) 



tan 6 = „^ I ZZTZTT^Jp > (44) 



For this step to be admissible we had to use the fact that cos 9, sin 9 ^ 0, which is true by assumption. 
Furthermore, we square and add (42a) and (42b), and (42c) and (42d), and divide the resulting equations. 
This yields the formula 

sin 2 + sin 2 \ tan 2 9 
cos 2 4> + cos 2 x tan 2 # 

which determines tan 5 up to a sign from known quantities <fi, x an d 6. However, we require [see (38)] 
that sin a, sin/3 always be nonnegative. Then the second equation in (43) implies that the sign of tan 6 
is the negative of the sign of tan x- Thus, 5 is determined up to multiples of n. Now, since sin/3 must be 
nonnegative, eq. (42b) determines the sign of cos<5, and therefore 5 is uniquely determined in [0, 2ir). 

Now that 5 is determined, a and [3 are uniquely determined by the system (42). In fact they are 
determined to lie in the interval [0, 7r), for, we have used the condition sin/3 > as well as the positivity 
of in (43) for the determination of S. 

Now let us examine case II.): Let 9 — 0. We first treat the case </>^ {0, f ,7r, ^} . Then the system 
(42) becomes 

cos 5 cos (3 = cos 4> , (45a) 
cos 5 sin/3 = , (45b) 
sin 5 cos a = sin (p , (45c) 
sin S sin a = . (45d) 

Since none of sm<j), cos</> is zero, eqs. (45a) and (45c) imply that sine), cosS ^ 0. But then eqs. (45b), 
(45d) imply that sin a = sin/3 = 0, and due to (38) it follows that a = (3 = 0. In this case, (45a) and 
(45c) imply that 5 = <fi. This gives the minimal decomposition (40a), which coincides with form (39a) for 
the choice of parameters (40b). 

Now we briefly discuss the cases where <j> <G {0, f ,7r, We proceed along the same lines as above, 
that is to say, using the system (45) as our starting point. For cos <fi — we arrive at the decomposition 

tf = R(-f)D(0)R(-§) , (46a) 

where a = and (3 is left unspecified. This indeterminacy corresponds to the non-uniqueness of the 
parametrisation (9,(f),x), eq. (22), when the parameter 9 takes values 0, § in the boundary of the 
parameter space. Eq. (46a) contains three matrices for (3^0, but only one for (3 — 0, since in this 
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case R(0) = 12- The choice (3 = therefore defines the minimal decomposition which is unique if (3 is 
restricted [eq. (38)] to lie in the admissible range [0,7r). This choice of parameters coincides with (40b), 
and the resulting form for U coincides with (40a). 

Analogously, for sin (f> = we have the decomposition 

C/ = R(-f)D(0)R(f) , (46b) 

with (3 = and a left unspecified. For general a, (46b) provides a decomposition in terms of three 
matrices. A minimal decomposition, in terms of one matrix only, is obtained for the choice a — 0, which 
is unique again, since a is restricted to lie in [0, 7r). Again, the parameter values a = (3 = corresponding 
to the minimal decomposition coincide with (40b), and the resulting minimal form is (40a). 

Now we examine case III.): Let = | but \ & {0; § i 71 "; %"} at nrs t. Then the system (42) becomes 

cos 5 cos (3 — , (47a) 

cos (5 sin (3 = cos x > (47b) 

sin S cos a — , (47c) 

sin S sin a = — sin \ ■ (47d) 

Since sin x, cos x ^ by assumption, we must have sin <5, cos 5 ^ from (47b, 47d). But then (47a, 47c) 
imply that cos a = cos (3 = 0, or a = (3 = |, due to (38). Then (47b, 47d) imply that <5 = — x, so all 
angles are uniquely determined, and the decomposition (41a) follows. 

Now we discuss the remaining cases, x € {0, f , 7r, 4^}. This proceeds exactly as for eqs. (46a, 46b). 
For cos x = we obtain the decomposition 

C/ = R( _£Z^ ) D(-x)R(^) , (48a) 

with a = \ and (3 left unspecified. This becomes minimal for (3 = ^, corresponding to (41b), and leads 
to eq. (41a). For sinx = we have 

^ = R( _o d ^ )D(x)R( o z ^ ) 5 (4gb) 

with (3 — f and unspecified a. The minimal decomposition is obtained for a = ^, corresponding to 
(41b), which leads to 

f/ = R(-f)D( X ) . (48c) 
However, in the case at hand we have D(x) = D(— x), so that the minimal form is again equal to (41a). 
This finishes the proof of theorem 4.1. ■ 



5 The Poincare Equivalence Theorem 

We are now in a position to deduce the Poincare Equivalence Theorem from our decomposition theorem 
4.1. Let us assume that U is the Jones matrix of a non-absorbing optical medium. As discussed before, U 
is unitary and may be assumed to be unimodular, possibly after discarding a global phase. Then U has 
a unique decomposition (39a) or (40a) or (41a). For the purposes of the Poincare Equivalence Theorem, 
we use the decomposition (39a) in terms of three matrices, so that 

E/ = R(-0D(J)R(C) • (49) 

Here we have introduced two angles 

£=— ' C = — ' (50) 

and it is understood that £ and £ are always chosen according to the minimal decompositions (39a, 40b, 
41b). In the next step we take into account that the diagonal matrix T>(6) is formally identical to the 
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Jones matrix of a linear retarder with fast axis along the x-direction and a phase lag of —28 of the slow 
with respect to the fast axis, 

D(5) = Jnn(0, -28) , (51) 
according to (10). We can now rewrite (49) as 

U = R(-0 J lin (0, -25) R(0 R(-£ + C) - Jiin(£, -25) R(C - £) (52a) 
= R(-e + C)R(-C) Jiin(0,-25)R(C) = R(C-O Jiin(C,-25) ■ (52b) 

Here we have used (15) and the group property (13) of the rotation matrices. 

In eqs. (52), the rotation matrices have been used to transform the Jones matrices of linear retarders 
at different angles of their fast axes into each other. However, in eq. (20) we pointed out that a rotation 
matrix R(— ^) is formally identical to the Jones matrix of a rotator J ro t(— 3r)- Then the right-hand sides 
of eqs. (52a, 52b) may be rewritten as 

U = J U n(^,-26) J rot (-/3) (53a) 
= J rot (-/3) Jh„(^,-25) . (53b) 



From the proof of theorem 4.1 we know that the decompositions (52, 53) are unique provided that 
the parameters (9, 4>, x) in the parametrisation (22) of U are restricted according to (29). 

If one of the parameter sets in (29) is taken, the matrices in (49) are no longer uniquely determined. 
From the proof of theorem 4.1 we know that in these cases either a or j3 remains unspecified. However, 
the phase retardation A = — 26 is always well-defined. It is easy to translate the condition that (9, 4>) 
or (6, x) take one of the values in (29a, 29b), respectively, into a condition for the optical parameters of 
the linear retarder and rotator in eqs. (53). This can be done with the help of (46a, 46b, 48a, 48b), or 
also (40a) and (41a). It is easy to check that the decomposition is non-unique if and only if the phase 
retardation A of the linear retarder in the equivalent optical model is an even or odd integer multiple of 
7r. Physically this means that the retarder in these cases can be realized by a stack of half-wave plates. 

In eqs. (52), (53) we have proven the 

Theorem 5.1 (Poincare Equivalence Theorem). For every non-absorbing passive optical medium 
(represented by a unimodular unitary matrix) there exists an optically equivalent model which is built 
from one linear retarder and one rotator only. These elements may be arranged in any order, and in both 
cases the phase lag A = — 25 is the same. 

If A is not a multiple of ir, a and (3 are uniquely determined. 

If A is a multiple of it ( corresponding to a stack of half-wave plates ), one of the angles a or (3 is left 
unspecified. However, a and [3 may still be chosen according to (40b) and (41b), corresponding to the 
minimal decompositions (40a) and (41a). 



6 Jones vectors, Stokes parameters, and the Poincare sphere 

In eqs. (2, 3) we have defined the polarization form and the associated Jones vector of a completely 
polarized beam of light. We mentioned that these quantities are defined only up to a global phase, so 
that Jones vectors v and e l * v are equivalent. 

An alternative description of the light beam in terms of the so-called Stokes parameters removes this 
ambiguity: The reason is that, in contrast to the Jones vector, the Stokes parameters are quadratic 
in the (complex) field strength and can, in principle, be determined by four intensity measurements, in 
conjunction with a linear polarizer and a quarter-wave plate, or equivalents [2, 31]. The Stokes parameters 
can be motivated by observing that, for a wave propagating in the z-direction, the scalar products 
d • E = ai and e 2 • E = a 2 are the complex amplitudes of electric field strength linearly polarized in 
the x- and y-direction, respectively. The phase difference between these components accounts for the (in 
general) elliptical polarization form of the light beam. Then the Stokes parameters arc defined in terms 
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of these quantities as 



So = |ei»E| + |e 2 • E| =a\ + a\ , 

I Till 2 I ttv|2 2 2 

si = |ei • E| — |e 2 • E| = a 1 — a 2 , 

s 2 = 2 %> [(ei • E)* (ei • E)] = 2 ai a 2 cos(5 2 - Si) , 

S3 = 2 3m [(ei • E)* (ei • E)] =2 ai a 2 sin(<5 2 - Si) . 



(54) 



The electric field vector of the polarized light beam, expressed by (3) or (54), describes an ellipse whose 
principal axes are not in the x- and y-direction in general. The first principal axis, with length 2a, makes 
an angle ip with the x-axis, while the second principal axis, being perpendicular to the first, has a length 
26. The shape and orientation of the ellipse can be described by the angle \ such that tan% = ±b/a. It 
is then a matter of elementary computation to show that the Stokes parameters (54) can be expressed in 
terms of these angles as [2] 

si = so cos 2% cos 2tp , 

S2 = so cos 2\ sin 2ip , > (55) 
s 3 = s sin 2\ ■ 

Evidently, the system (55) can be interpreted as the Cartesian coordinates of a point on a sphere of 
radius so, the so-called Poincare sphere [27]. The significance of the Poincare sphere lies in the fact that 
its points are in 1-1 correspondence with the distinct states of polarization of a completely polarized 
light beam with constant intensity; the ambiguity in the Jones matrix formalism, in which a continuous 
infinity e'* v of Jones vectors corresponded to physically equivalent polarization forms, has been removed. 
One can show that right- (left-) handed circular polarization is represented by the north (south) pole on 
the Poincare sphere; whilst the states of linear polarization correspond to the points in the equatorial 
plane. All other points describe a general elliptical polarization. Each of the equatorial points represents 
a different angle of the plane of linear polarization; as a consequence, an optical rotator, changing just 
this angle, but preserving the linear polarization form, may be regarded as a rotation i?3 about the z- 
axis on the Poincare sphere. Similarly, an optical retarder, capable of transforming a linear polarization 
state into a general elliptical one, will correspond to a rotation about an angle 2-0 whose axis lies in the 
xy-plane of the coordinate system. However, any such rotation can be obtained as a sequence i? 3 i? 2 R% 1 
of rotations about the z-axis, the momentary y-axis, and the new z-axis again. 

More generally, it can be shown that the action of any loss-less optical element corresponds to a proper 
rotation on the Poincare sphere. From this fact, the Poincare Equivalence Theorem can be inferred 
immediately: For, let U denote the unitary matrix of the optical clement in the Jones matrix formalism, 
and let R(U) denote the corresponding rotation on the Poincare sphere. We know that every proper 
rotation can be decomposed into three Ruler rotations, which may be chosen as a sequence R' 3 i? 2 R3 of 
three rotations about the momentary z- and y-axes. From the previous paragraph we know that the 
factors in this decomposition represent the actions of rotators and linear retarders, so that the sequence 
R' 3 R2 i?3 corresponds, in the Jones matrix formalism, to a sequence of optical elements R' D R, which is 
of course just eq. (39a). In this way, the decomposition of 5?7(2)-matrices into 'basic building blocks' as 
in the Poincare Equivalence Theorem can be traced back to the Euler decomposition of proper rotations 
on the Poincare sphere. 

The parameters in (54) are classical quantities, of course. It is possible to define associated quantum 
Stokes parameters together with a quantum Poincare sphere; these concepts are reviewed in [16], for 
example. In [15], an application of the quantum Stokes parameters to squeezed light has been presented. 



7 Relation between Characteristic Parameters and the Poincare- 
Equivalent Model 

In section 3 we have shown that every optical device which (in the framework of Jones calculus) can be 
represented by a unitary unimodular (2,2) matrix U possesses two orthogonal directions, wi and w 2 , 
such that linearly polarized light whose plane of polarization on entry into the apparatus is Wi or w 2 
will emerge linearly polarized again, with polarization directions w[ and w' 2 - Furthermore, we saw that 



13 



linearly polarized light entering at the two primary characteristic directions emerges with different phases; 
an expression for the relative phase difference was given in (36), but we were not able to determine it 
uniquely. 

On the other hand, in theorem 5.1 we showed that every such device was optically equivalent to a 
succession of linear retarders and rotators. Since a linear rctarder has two preferred spatial directions (the 
fast and slow axis) and a well-defined phase retardation between them, this suggests that the characteristic 
directions might be closely related to the fast and slow axes, while the phase difference (36) might be 
identified with the retardation of the retarder. 

We now show that this is indeed so: We first treat the unique cases, i.e. parameters (6, <fi, \) satisfy 
(29), or equivalcntly, the phase lag A = — 25 of the equivalent linear retarder is not a multiple of tt. In 
this case the linear retarder has precisely two preferred directions, namely the fast and slow optical axis. 
Only at these directions preserves the linear retarder the state of linear polarization; any other direction 
of linear polarization on entry will emerge as elliptical polarization form. On using the decomposition 
(49) it is then clear that the primary characteristic directions are determined by the condition 



R(()w m = e m , TO =1,2 , 
ei = (l,0) T , e 2 = (0,l) T . 



(56) 



Therefore we must have 

toi = R(-C) ei 



sin (y ) 



, w 2 = R(-C) e2 



Mj ' 



cos I 



(57) 



where we have defined that w\ be the direction corresponding to the faster phase velocity. Similarly, the 
secondary characteristic directions must be 



4 = R(-0ei 



cos 
sin 



> w 2 = R(-C) e 2 



in(^) 



-sm 
cos 



(58) 



Again, we have defined that w[ be the secondary characteristic direction corresponding to the faster phase 
velocity. A comparison between (57) and (23), and (58) and (24) now shows that we should set 



a-B 

71 = -a* = 



C , 72=7i + f , 



and 



where we have used (50). 



7i = ^ = e 



7 2 = 7i + | 



2 ' 



(59a) 
(59b) 



Light which travels through the medium with linear polarization along w\, w 2 acquires a phase 
(relative to a global phase <J> g ) of 

$ x = 5 , $ 2 = -5 . (60) 



Then the total phase difference between directions u>i and w 2 is 

A = $ 2 - $i = -25 . 



(61) 



The arguments leading to eqs. (57-61) clearly are based on intuition rather than rigorous derivation. 
In order to be rigorous we have to show that these equations are consistent with the results of theorems 
3.1, 3.2 and eq. (36) in section 3. More precisely, we must show that the angles 7 m as defined in (59a) 
satisfy the equations (25), while "f' m as defined in (59b) must satisfy (27). This can be done easily, on 
using the assumption that parameters (9,<p,x) satisfy (29): To show that the first statement is true it is 
sufficient to expand tan(a — (3) in terms of sines, cosines of a and (3, and multiply both nominator and 
denumerator of the resulting fraction with cos 5 sin 5; it is at this point where we need the condition that 
cos 5, sin 5 ^= 0. The result can be easily recast in the form of eq. (25). The second statement linking 
(59b) to (27) is confirmed in the same way. Finally, eqs. (33a, 33b) can be derived immediately from the 
system (42). We note that we have removed the indeterminacy concerning the assignment between $ m 
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and ±S by decreeing that the subscript 1 should refer to the fast axis of the equivalent retarder. This 
assignment is now well-justified since it is based on the distinction between different values of a physically 
measurable quantity, namely the phase velocity along the axes. Previously, in section 3, we did not yet 
possess the physical justification for such a deliberate choice. 

We now turn to discuss the non-unique decompositions: From section 5 we know that in these cases 
the phase retardation of the equivalent linear retarder is always a multiple of 7r; in particular, it is always 
well-defined. What is left unspecified is one of the angles a, j3. This is physically reasonable, for, whenever 
A = kit, k G Z, the retarder transforms any linearly polarized light into linearly polarized light again, 
and hence any direction (before entry into the medium) is a primary characteristic direction (and any 
direction behind the medium is a secondary one). This is clearly consistent with the eqs. (25, 27): In 
the comment following (29) we remarked that (25, 27) remained ill-defined if parameter values (29) were 
taken. This simply means that there are no preferred directions in this case, and any two orthogonal 
directions before and behind the medium may serve as primary and secondary characteristic directions. 

We finish this section by writing down the decompositions (53) of the equivalent optical model in 
terms of the characteristic parameters: In the case of unique decompositions, eqs. (49) and (53) read 

U = J ro t(-7i) -WO, A) J rot (7i) (62a) 
= Jiin(7i,A) J ro t(7i-7i) (62b) 
= Jrot(7i-7i) Jiin(7i,A) , (62c) 

where we have used (59) and (61). For 6 = 0, the minimal decomposition (40a) can be expressed as 

U = Ji in (0, A) , A = kit , k e Z . (62d) 

For 9 = | , the minimal decomposition (41a) can be expressed as 

U = Jrot(-f) Jlin(A) , A = kTT , fceZ . (62e) 



8 Basic Ideas of Integrated Photoelasticity 

In this section we want to outline, in a single paragraph, how optical media with spatially dependent 
dielectric tensors naturally emerge in Photoelasticity. 

Consider a medium (such as glass or certain polymers) which is optically homogeneous and isotropic, 
with a dielectric tensor e5ij when strain-free, but which becomes inhomogeneous and anisotropic when 
subject to an external load. The dielectric tensor then is related to the stress tensor ct^ of the material 
according to the stress-optical law 

e ij = 6 % + Ci a ij + ^2 tr(er) 5 t j , i, j = 1, 2, 3 , (63) 

where C\, C2 are stress-optical constants. The medium is assumed to occupy a bounded region in space. 
Polarized light is sent through the specimen at many different angles, and from the change of the state of 
polarization the characteristic parameters of the model, for a given direction of the incident light beam, 
are retrieved. Light is assumed to propagate through the medium along straight lines, an approximation 
which is acceptable as long as the stress-induced birefringence is weak. Under these circumstances, the 
Maxwell equations describing the propagation of a plane wave with locally varying state of polarization 
along the z-direction of the medium can be cast into the form 

dz\v 2 ) 2c^i\e2i ^22 - e J \v 2 J ' 

where v\(z), V2(z) are the components of the spatially dependent Jones vector for fixed x and y, and 
tij(z), i,j = 1,2 are the components of the projection of the local dielectric tensor onto the xy-plane, at 
position (x,y, z). The point (x,y) in the xy-plane characterizes the location of entry of the light beam 
(which is assumed to have very small diameter) into the medium. Eqs. (64) may be expressed in an 
equivalent way by introducing a unitary propagation matrix U(z) such that 

v(z) = U(z,Zi)v(zi) , (65) 
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where Zi is any point before the medium. This matrix must satisfy a differential equation analogous to 
(64), 

■%-U(z,Zi) = iH(z)U(z,Zi) , (66) 

where the matrix H is given on the right-hand side of (64). If Zf is a point behind the medium, the matrix 
U(z/,Zi) is obtained by integrating (66); it determines the change of polarization form of the light beam 
passing through the medium between points (x,y, Zj) and (x,y,Zf) along a straight line. Measurements 
of this change in the polarization determine the characteristic parameters, and in turn, the equivalent 
model in either of the forms (62), for the given ray of light. Repeating the procedure for many rays at 
different angles yields data sets which contain the characteristic parameters for any possible direction of 
the light beam. From these data we attempt to reconstruct the local dielectric tensor in the interior of 
the medium. This problem, called three-dimensional Photoelasticity, is still unsolved in its most general 
form, and the Poincare Equivalence Theorem as given in one of the forms above may play an important 
role in this goal. 



9 Summary 

In this paper we have proven a decomposition theorem for two-dimensional unitary unimodular matrices 
which represent a non-absorbing passive optical medium in the framework of the Jones calculus. From 
this decomposition theorem, the Poincare Equivalence Theorem can be derived as a corollary. The latter 
states that every non-absorbing optical medium is optically equivalent to a succession of linear retarders 
and rotators, in any order. Except for few cases, the equivalent linear retarder and rotator are uniquely 
determined from data characterizing the original medium. We have shown that the parameters of these 
equivalent elements have a simple interpretation in terms of measurable quantities, the so-called Char- 
acteristic Parameters of an optical model, which have been known in the engineering literature for some 
time. The relation between these Characteristic Parameters and the Poincare-equivalent optical model 
has been fully clarified. Finally, we have illustrated the context in which the Poincare-equivalent model is 
expected to become a useful tool: Within the framework of three-dimensional Photoelasticity, the decom- 
position theorems proven in this paper may turn out to be important to the goal of reconstructing local 
dielectric tensors from measurements of the global change of the polarization state of light transmitted 
through an inhomogeneous and anisotropic medium with spatially varying refractive tensor. 
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